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Abstract 

We relate ordinary and skewed parton distributions to soft overlap contributions 
to elastic form factors and large angle Compton scattering, using light-cone wave 
functions in a Fock state expansion of the nucleon. With a simple ansatz for the 
wave functions of the three lowest Fock states we achieve a good description of 
unpolarised and polarised parton distributions at large x, and of the data for the 
Dirac form factor and for Compton scattering, both of which can be saturated with 
soft contributions only. Large angle Compton scattering appears as a good case 
to investigate the relative importance of soft and hard contributions in exclusive 
processes which are sensitive to the end point regions of the nucleon wave function. 
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1 Introduction 



The recent theoretical developments for real and virtual Compton scattering, which have 
lead to the introduction of skewed parton distributions^ (SPDs) |T], ^ |3|, have renewed 
the interest in the interplay between hard inclusive and exclusive reactions. In the light- 
cone approach the link between these classes of reactions is mediated by light-cone wave 
functions (LCWFs). Although this connection has been known for quite some time [|], [5| 
it has not yet been exploited practically. 

An important question in this context is the size of perturbative QCD contributions to 
exclusive reactions. There is general agreement that the conventional hard scattering ap- 
proach (see [[| and references therein), in which the collinear approximation is used, gives 
the correct description of electromagnetic form factors and perhaps other exclusive pro- 
cesses in the limit of asymptotically large momentum transfer. The onset of that asymptotic 
behaviour is however subject to controversy. It has turned out that for the electromagnetic 
form factors of the pion and the nucleon or for Compton scattering agreement between data 
and the perturbative contributions is only obtained if distribution amplitudes are employed 
that are strongly concentrated in the end point regions, where one of the parton momentum 
fractions tends to zero. Such distribution amplitudes have been proposed by Chernyak et 
al. |J on the basis of QCD sum rules, but their derivation has been severely criticised, cf. 
for instance J7|. At least for form factors but likely also for Compton scattering they lead 
to perturbative contributions which are dominated by contributions from the end point 
regions where the use of perturbative QCD is not justified M. In the case of the pion dis- 
tribution amplitudes concentrated in the end point region are now excluded by the CLEO 
data H on the 7r7 transition form factor, where they lead to perturbative contributions 
much too large in comparison with experiment Jl0|, [IT[] . In the case of the nucleon form 
factor it has been shown in Ref. Jl2|] that the inclusion of transverse momentum effects 
as well as Sudakov suppressions ||13[ in the perturbative analysis leads to a substantial 
reduction of the perturbative contribution which then is much smaller than experiment. 

There is another difficulty with distribution amplitudes concentrated in the end point 
regions: if they are combined with a plausible Gaussian transverse momentum dependence 
in a wave function and if from that LCWF the soft overlap contribution [II to the nucleon 
form factor is evaluated one obtains a result that exceeds the form factor data dramat- 
ically M. Such wave functions also lead to valence quark distributions that are much 



larger at large x than those derived from deeply inelastic lepton- nucleon scattering |p~5 |. 
Starting from all these observations and from the assumption of soft physics dominance, 
the authors of Ref. |TB[ derived a LCWF for the nucleon's valence Fock state by fitting 
its free parameters to the valence quark distribution functions and the form factors in the 
momentum transfer region from about 5 to 30 GeV 2 . The LCWF obtained in [TIJ is close 
to the asymptotic form and very different from the end point concentrated ones. Recently 
Radyushkin |TTJ generalised the overlap approach to large angle Compton scattering and 

lr rhe name skewed parton distributions has been proposed to amalgamate the different terms (nonfor- 
ward, off-forward, nondiagonal, off-diagonal) used in the literature for closely related quantities. 
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showed that soft physics, evaluated from LCWFs similar to the one used in [TBI], can account 



for high energy Compton scattering in the experimentally accessible kinematical region as 
well. It goes without saying that the soft contributions to form factors and Compton scat- 
tering are suppressed by inverse powers of the hard scales compared with the perturbative 
contributions, which will always dominate at very large energy and momentum transfer. 

The purpose of the present paper is firstly to extend the analysis of |16[ to higher Fock 
states in order to explore their importance relative to the lowest one, and secondly to 
include Compton scattering in the analysis, following Radyushkin's work [IT]]. In Sect. |2| 
we present some kinematics of the elastic form factor and of Compton scattering. We then 
give a general discussion concerning soft contributions and the essential conditions for a 
representation of form factors and other processes as an overlap of LCWFs (Sect. |]). Soft 
contributions to real and virtual Compton scattering arising form the handbag diagrams 
will be discussed in Sect. |J In the next section, Sect. ||, we introduce our parametrisations 
of LCWFs for the lowest Fock states. In Sects. |5], [Tj, § we respectively evaluate parton dis- 
tributions, form factors and large angle Compton scattering. As an extension of evaluating 
parton distributions in the Fock state approach we also calculate skewed parton distribu- 
tions (Sect. |). Since our LCWFs describe quite well the quantities mentioned before, our 
results for the skewed distributions may convey an impression how these functions look 
like. The paper ends with our summary (Sect. [Top . 



2 Kinematics 

To begin we give our notation for the elastic form factor and for Compton scattering and 
introduce several reference frames we will need later. 

2.1 The elastic form factor and Compton scattering 

The external momenta of the one- and two-photon processes 7*p — > p and 7*p — > 7p are 
denoted as shown in Fig. p] (a) and (b). We use the Mandelstam variables s = (p + q) 2 , 
t = A 2 , u = (p — q') 2 , and write Q 2 = —q 2 for the incoming photon virtuality in Compton 
scattering and m 2 = p 2 = p' 2 for the squared proton mass. Note that we write A (and 
not q) for the momentum transfer to the proton in the elastic form factor, reserving q for 
the incoming photon in the Compton process; this will be useful to display the similarities 
of the one- and two-photon processes. We denote the momenta of the active partons, i.e. 
those that couple to the photons by k and k', and for the parton-photon subprocess in 
Fig. [1] (b) we use Mandelstam variables s = (k + q) 2 , i = t and u = (k — q') 2 . Whenever it 
is necessary to distinguish the momenta of active and spectator partons we will label the 
active one with an index j and the spectators with an index i (i ^ j); outgoing momenta 
will always be indicated by a prime. 

In the various reference frames described below we introduce light cone variables v ± = 
(v°±v 3 ) / v2 and the transverse part vj_ = (v 1 , v 2 ) for any four-vector v and use component 
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Figure 1: Overlap diagrams for (a) the elastic form factor and (b) Compton scattering. 
Lines p and p' denote protons, k and k! quarks or antiquarks, and the horizontal lines 
represent any number of spectator partons. The small blob attached to the photon lines 
stands for the pointlike quark-photon coupling in (a) and for the two diagrams of Fig. |2| 
in (b). 



notation v = [v + ,v , vj_]. We finally define the ratios 
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(1) 



of plus-components; positivity of the energy of the final state proton and photon implies 
C < 1 and K) > 0. 

Let us take a closer look at the physical region of the variables t and (. In any reference 
frame we can write 
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using our definition of ( and the on-shell conditions for the proton momenta. With 
and ( < 1 we have 
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which imposes a minimum value on —t at given (. Note that this is independent of the 
process considered. 



2.2 A symmetric frame 

For reasons that will become apparent in Sect. ^ frames where A + = 0, i.e. ( = 0, play a 
special role in the context of overlap contributions. We shall use a frame where 
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which treats the transverse momenta of incoming and outgoing hadron in a symmetric way 
and presents the further simplification that A~ = 0. Note that t = — A.\ here. Condition 
(||) fixes the frame up to a boost along the 3-axis. For the elastic form factor one may 
take any frame satisfying (|J); in the case of Compton scattering a symmetric choice is to 
further impose p 3 + q 3 = 0. Note that for real Compton scattering this is just the cm. 
frame with the 3-axis along p + p', while with a virtual initial photon it does not coincide 
with the centre of mass. For the photon momenta we have 
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(t + Q 
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2t 
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2t 



with 
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t + Q 2 
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s + u — 2m 



(5) 



(6) 



5 — 11 + 2^1 s(u — Mi) (t — t)/t 

where t = t corresponds to forward and u = U\ to backward scattering in the photon- 
proton cm. We shall in the following refer to this frame as the "symmetric frame" . 



2.3 The photon-proton cm. 



As we will see in Sect. |3.1.2| the symmetric frame just described is not suitable for our 
discussion of deeply virtual Compton scattering (DVCS). In that case we will use the cm. 
frame with the 3-axis pointing in the incoming proton direction, where 
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and again p 3 + q 3 = 0. The non-vanishing plus component of the momentum transfer is 
characterised by the skewedness parameter (; the total momentum transfer to the proton 
reads 

+ (m 2 + A 2 



and its square is 
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2p+(l-C) 
Cm 2 + A 2 ± 



(9) 



(10) 



Notice that the relation (p^0[ ) follows from (0) alone and thus holds in any frame where 
p^ = 0. In the photon-proton cm. frame we have 

Q 2 -t(l 
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7] = ( — X 
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where 

X N = . ^ XBj , XBj = — — (12) 

l + v /l + 4x2 j mVg 2 2 P-V 

respectively denote Nachtmann's and Bjorken's variable. In the kinematical region of 
DVCS, i.e. when —t is small, and Q 2 and s are large (|TT| ) simplifies to £ ~ a;^ and 77 ~ 0. 



2.4 Frames for the hadron LCWFs 

The arguments of LCWFs are given as the plus-momentum fractions X{ and the trans- 
verse parts kjj of parton momenta in a frame where the transverse momentum of the 
corresponding hadron is zero. We will call those systems "hadron frames" and refer to 
transverse parton momenta in an appropriate hadron frame as "intrinsic" transverse mo- 
menta. 

The transformation from a given frame to a hadron frame can be achieved by a "trans- 
verse boost" (cf. e.g. [ 1 8 1 ) which leaves the plus component of any momentum vector a 



unchanged, and which involves a parameter b + and a transverse vector bj 
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Starting for instance from the symmetric frame of Sect. the choice b + = p + , b 
— A^/2 transforms the momenta of the incoming hadron and its partons as 
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(14) 



where we suppressed the minus components of the parton momenta, whose expression we 
will not need. This is an appropriate frame to read off the arguments of the LCWF of the 
incoming hadron as Xi and kj_j = kj_j + Xi Aj_/2. The analogous boost with the choice 
b + = p + , bj_ = +A.±/2 relates the symmetric frame with a frame appropriate for identifying 
the arguments of the LCWF for the scattered hadron as x\ and k' ±i = k' ±i — x\ Aj_/2. 

Incoming and outgoing parton momenta in the overlap contributions Fig. [I] are related 
by k\ — k i (i 7^ j) for the spectator partons and kj = kj + A for the active parton which 
takes the momentum transfer in the scattering. Using the transformations between the 
symmetric frame and the in/out-hadron frames established above we can directly express 
the LCWF arguments for the outgoing hadron (denoted by a hat) in terms of the ones for 
the incoming hadron (denoted by a tilde): 



x i 



x. 



Xa 



x, A 



for i ^ j , 



Kj + (1 " xj) A_ 



(15) 



where we could have dropped the hat/tilde notation for the momentum fractions which 
are not changed by the boost (ITT 
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For the case of deeply virtual Compton scattering the photon-proton cm. frame intro- 
duced in Sect. |2.3| is already the appropriate hadron frame to identify the arguments of the 
LCWF of the incoming proton. By the boost ([13]) with the parameter values h + = (1— C) P + , 
b± = A.± one obtains the momenta in the corresponding frame where the outgoing proton 
has zero transverse momentum. LCWF arguments for the outgoing proton (denoted by a 
breve) are related to the ones of the LCWF of the incoming proton as 

k^ = k ±i + i^|A ± , (16) 

where according to its definition the plus momentum fraction in the LCWF of the scattered 
proton is taken with respect to p' + = (1 — () p + and not to p + . We notice that for ( = 
Eq. ( |T6"D takes the same form as flT5|). 



X ; 



i-C 

Xj -C 



c 



3 The theory of soft contributions 

In this section we are concerned with soft overlap contributions to hard exclusive processes. 
They are contributions where only some of the partons in the external hadrons are active, 
i.e. participate in a hard scattering, while the other partons remain spectators. 



3.1 Bethe-Salpeter and light cone wave functions 

The evaluation of overlap contributions in terms of light cone wave functions requires 
some care. An example is the Drell-Yan overlap formula [14] of the elastic form factor, 
for which Isgur and Llewellyn Smith M observed that different results are obtained in 
different reference frames. Sawicki |19| has shown the origin of this discrepancy: in certain 
reference frames there are overlap contributions which are not contained in the Drell-Yan 
formula; when they are taken into account Lorentz invariance is restored. We shall first 
review Sawicki 's arguments J 



20] for the form factor and then investigate the case of 



Compton scattering. 



3.1.1 The elastic form factor 

Our starting point to obtain the overlap formula for the form factor is the diagram of Fig. [I] 
(a) in the framework of equal-time quantisation and covariant perturbation theory. The 
hadron-parton vertices, represented by the large blobs in the diagram, are described by 
Bethe-Salpeter wave functions ^bs- For simplicity we consider a scalar hadron coupling to 
two scalar partons, so that there is only one spectator line in the diagram Fig. [1] (a). We 
further work in a toy theory where the hadron has a pointlike coupling to the two partons; 
to leading order in the coupling constant the wave function ^Bs(k) of the hadron with 
momentum p is then given by the coupling times the free propagators for the partons with 
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momenta k and p — k. In general (and in particular for QCD) ^B,s{k) will have a more 
complicated analytic structure in the virtualities k 2 and (p — k) 2 involving branch cuts in 
these variables. Their discussion is beyond the scope of this paper and we only retain the 
propagator poles in these variables. This will be sufficient to exhibit the points we want 
to make. 

The aim is now to perform the loop integration over k~ in Fig. [I] (a) so as to reduce 
^Bs{k) and ^bs{^') to LCWFs. For this we use that up to a normalisation factor a 
LCWF is obtained from the corresponding Bethe-Salpeter wave function, say by 
the integral J dk~ ^Bs{k) at fixed k + and kj_. Note that this relation does not only hold 
in frames where the hadron has zero transverse momentum; with (|H|) we see in particular 
that this integral is invariant under transverse boosts. The ^"-integration in Fig. [I] (a) is 
readily performed using Cauchy's theorem since the analytic structure of the diagram is 
given by the propagator poles in our model. Writing 



and using that the poles in k 2 , (k — p) 2 , (k + A) 2 are situated just below the real axes 
in these variables we see that the propagator poles are below or above the real k~-axis, 
depending on the value of x. For definiteness we now take ( > 0, where we have the 
following cases: 

1. For x > 1 and for x < all poles are on same side. Closing the integration contour 
in the half plane where there are no singularities one obtains a zero integral. 

2. For 1 > x > ( we pick up the pole in (p — k) 2 alone when closing the integration 
contour in the upper half plane. The diagram is then given by the propagators of k 
and k + A, evaluated at the value of k~ where p — k is on shell. Applying an analogous 
argument to the integral J dk~ ^Bs{k) we find that a LCWF can be written as the 
hadron-parton coupling times one parton propagator, evaluated at the value of k~ 
where the other propagator is on shell. As a by-product one finds that the plus 
momentum fractions of the partons w.r.t. the hadron are always between and 1, 
otherwise the integral is zero. In total we find that the diagram for the form factor 
is given by the product of two hadron LCWFs, as stated in the Drell-Yan formula. 

3. For ( > i > we can pick up the residue at the pole in k 2 , or alternatively the sum 
of residues for the poles in (p — k) 2 and (k + A) 2 . In the term where k (or k + A) is on 
shell both partons in the hadron p' (or p) are both off-shell, which cannot be rewritten 
in terms of a LCWF. Note that for ( > x > the parton that has been struck by the 
photon has negative plus-momentum fraction x — (, which does not correspond to a 
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(17) 
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parton going into hadron p'; a situation that clearly cannot be expressed through a 
LCWF. This contribution is missing if one naively writes down the Drell-Yan formula 
in a frame where ( 7^ 0: here is the origin of the paradox observed in 

We thus see that in order to obtain an overlap representation in terms of two-parton 
LCWFs for each hadron we need to go to a reference frame where (, or in other words, A + 
is zero. Such a frame was in fact chosen in the original work by Drell and Yan.f] 

The argument goes along the same lines if one has more than one spectator and takes 
a Bethe-Salpeter wave function with only the propagator pole in each parton. Let us label 
the active parton (i.e. the one hit by the photon) k\ and the spectators k 2 , &n-i, 
k n = p — ki — k 2 . . . fc n -i and first perform the integrations over k 2 , . . . , k~_ x to put n — 2 
spectators on shell while not doing anything to the active parton. Then we are left with 
one active parton and a cluster of spectators, and the situation for the integration over k^ 
is as above. 



3.1.2 Compton scattering 

We shall now see that the two-photon process of Fig. |1] (b) involves a new difficulty Let 
us take again our toy model of a hadron with a pointlike coupling to two partons. To 
leading order in the electromagnetic coupling the parton-photon vertex is given by the two 
diagrams of Fig. |2|. Compared to the form factor case we thus have an extra propagator in 
the overall process, corresponding to a squared momentum s = (k + q) 2 or u = (k — q') 2 , 
so that fllTD is completed by 



{k + q) 2 + (k_L + q ± 



,2 



in the s-channel and 



k~ = - — / V -q~ (18) 

2p+(x- V ) +q 1 yj 

in the w-channel diagram. 

In the case 1 > 77 > ( > one has the possibilities listed in Tab. |1] to pick up propagator 
poles in the fc~-plane. Proceeding in the same way as in the form factor case we see that 
only in the region 1 > x > 77 we obtain an expression in terms of the two-parton LCWFs 
for both hadrons. In all other regions we have further contributions, either from a parton 
attached to one hadron but not the other (k or k + A) or not attached to a hadron at all 
(k + q or k — q'). 

2 In a recent paper pi] this contribution has been rewritten in terms of the LCWF for the hadron p 
containing partons with momenta k and — (k + A) plus the hadron with momentum p', and the (trivial) 
LCWF to find the hadron with momentum p' in the hadron p' . It would be interesting to explore this idea 
in the context of Compton scattering, but this shall not be done here. 

3 In a frame with £ = there can be finite contributions of the type discussed in point 3 if the integrand 
in the interval £ > x > becomes singular for ( — 0. This happens for the minus component of the parton 
current ]2l|, [2^], which we shall not use in our applications; the plus component of the current does not 
exhibit this phenomenon. 
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H 



q' = q - A 



k' = k + A 



X 



q' — q — A 



(b) 



fc' = fc + A 



Figure 2: (a) s-channel and (b) w-channel diagram for quark-photon or antiquark-photon 
scattering. 



diagram 


region 


propagator pole in 


s- channel 


1 > x > C 
C > x > 

> X > ( — T] 


p — k or k, k + A, k + q 
k, k + q or p — k, k + A 
k + q or p — k, k + A, k 


u- channel 


1 > X > T] 
1] > X > ( 

C > x > 


p — or k, k + A, k ~ q' 

p — k, k — q' or k, k + A 

k or p — k, k — q' , k + A 



Table 1: Possibilities to pick up propagator poles in the fc -integration for the case 1 > 
rj > C > 0. 

The situation is analogous in other cases than 1 > rj > ( > 0, and also if there is 
more than one spectator. Notice that in general we cannot find a frame where r] = ( = 
to solve our problem: if 77 = then q' 2 = implies q' ± = 0, and if also ( = then 
t + Q 2 = -2q ■ q' = which is a special kinematical situation. 

At this point we look beyond our toy model and remember that we want to evaluate soft 
overlap contributions in QCD, which involve the soft parts of the hadron wave functions, 
not the hard parts that are generated perturbatively ||. We will now see that in certain 
cases we can obtain an approximate expression for the soft overlap contribution that in- 
volves only the LCWFs of the two hadrons. To this end we first chose a frame with £ = 
so as to eliminate the interval ( > x > 0, as we did for the form factor. It turns out that 
with appropriate external kinematics the contributions from the poles in s and u go with a 
highly virtual parton in at least one of the two hadrons. Since large parton virtualities are 
strongly suppressed in the soft parts of the hadron wave functions (they constitute their 
hard parts) we can neglect these pole contributions, restrict x to the interval from to 1 
and only take into account the contribution from the pole in (p — k) 2 , which just leads to 
an expression with two hadron LCWFs as in the form factor case. 

To see when this is the case we write 



s + Q 2 



x(s + Q 2 



m 



+ k 2 



10 



- 2(k_L - xp±) ■ q_L - \x 2 m 2 - k 2 - (kj_ - xp±) ■ (k_L + xp^)} , 

x L > 

= x{u — m 2 ) + k 2 

+ 2(k ± -xp ± )-q' ± + ^{x 2 m 2 -k 2 -(k ± -xp ± )-(k ± + xp ± )} (20) 



and make the hypothesis that the soft hadron wave functions are dominated by intrinsic 
transverse parton momenta kj_j satisfying k 2 Li /x i ^, A 2 (this is implemented in our ansatz 
for the LCWFs in Sect. fD, where A is a hadronic scale in the GeV region, and by parton 
virtualities in the range \k} \ ^A 2 . From now on we concentrate on two cases. 

Large angle Compton scattering (large s, —t and —u) 



We now work in the symmetric frame of Sect. |2.2| . Let us for a moment stick with the 
case where there is only one spectator parton; the expressions k 2 Li /x i and k'fjx^ for this 
spectator in the initial and final state hadron (cf. Sect. |2.4| ) can then be rewritten in terms 
of the active parton momenta k and k! . For their sum we obtain 

(k ± - ^ + (k- - *pl)* + 2(k, + Ax/2? < A2 



1 — x 1 - ./• " 1 — X 

which implies 



|1 -x|<A 2 /(-t) , (kx-xpxI^A 2 /^. (22) 

In the case of several spectators the argument goes along the same lines, now summing 

kii/^i + ^h/x'i over an spectators. 

We remark in passing that a restriction to intrinsic transverse momenta k 2 u < A 2 instead 
of kjj/a^^A 2 would not be enough to ensure small parton virtualities in the hadrons: 
instead of ( Ell ) we would then only have (kj_ — xpj_) 2 + (k' ± — xp' ± ) 2 ^A 2 , which gives 
1 1 — a; | ^ A/a/— t and |kj_ — xp^| ^5 A, and in particular |kj_ + Aj_/2| ^ A. From fc 2 — k' 2 = 
2Aj_ • (kj_ + Aj_/2) we see that then at least one of the parton virtualities would be of 
order A^/— t and not A 2 . 

With (|2~0"|), (p^) and (||), @ we have s ~ s and u ~ m up to corrections of order 
A 2 (t ± Q 2 )/t, provided that both s and — u are large on a hadronic scale.0 This implies 
that in order for s or u to have a pole at least one parton must have a large virtuality or 
intrinsic transverse momentum, so that following our above remarks we can neglect these 
pole contributions. Note that apart from —t one also needs — u large: when the latter 
becomes too small the propagator u can easily become soft and it is no longer justified to 
neglect its pole contribution (which one may relate to the soft, hadronic part of the final 
state photon). Similarly one can see that s must be large, too. 

The physical situation clearly is that of a hard photon-parton scattering and the soft 
emission and reabsorption of a parton by the hadron, similar to the familiar handbag 

4 One may admit a two-scale regime A 2 ^ — t <C Q 2 provided that s and — u are also of order Q 2 . 
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diagram for inclusive deeply inelastic scattering (DIS) or for DVCS.f] In the hard scattering 
one can approximate the parton momenta k, k' as being on shell, collinear with their parent 
hadrons and with light cone fractions x — 1. This also provides another point of view on 
neglecting the s and u pole contributions: approximating k~~ with the value for which the 
partons are on shell in the hard scattering we have a ^"-integral where only the parton 
lines directly attached to hadrons provide a ^"-dependence. This is just as in the case of 
the elastic form factor, and thus we have the same situation for expressing the amplitude 
in terms of LCWFs as described in Sect. |3 . 1 . 1| . 

At this point we can also understand why in the conventional hard scattering mecha- 
nism [f|] (and also in the modified one of Botts, Li and Sterman |13|) one always obtains an 
expression involving hadron LCWFs, irrespective of the reference frame used. The reason 
is that in the corresponding diagrams the parton lines from each hadron are directly at- 
tached to a hard scattering subprocess, where the minus components of their momenta can 
be approximated with their values for which the partons are on shell. The corresponding 
^"-integration then only concerns the hadron-parton vertex alone and leads to a LCWF. In 
the case of soft overlap diagrams the situation becomes more complicated because spectator 
parton lines are "shared" by different hadrons, without undergoing a hard scattering. 



Deeply virtual Compton scattering (small —t, large Q 2 and s) 



In the kinematical region of DVCS, where — t ~ A 2 , we can no longer infer from (|2l|) that 
x must be close to one. Furthermore the factors (t ± Q 2 )/(2t) in (g), (|) are large, and 
the terms involving q ± , q' ± and 77 in (|20|) can thus be of order of the large scale,Qso that s 
or u may be zero even if the partons are near shell and nearly collinear with their parent 
hadrons. Our previous argument to neglect the pole terms in s and u then no longer works 
in the frame we have considered so far. There exist other frames with ( = 0, but one can 
show that 1] cannot be smaller than in (^) by solving the minimisation problem for 7] with 
an arbitrary axis defining plus components under the constraint £ = 0. 

We know however from the factorisation theorem of DVCS |3|, |23| that in a frame such 
as the cm. where the incident and the scattered hadron move fast to the right (and where 
C 7^ 0), the process factorises into a skewed parton distribution describing the soft coupling 
between partons and hadrons, and a hard photon-parton scattering calculated with the 
minus- and transverse components of k and k! replaced with zero. This factorisation is not 
realised in our symmetric frame of Sect. [2.2| , where the hadron momenta become slow of 
order y/—t in the DVCS limit. 

Using this factorisation in the cm. we can again neglect the pole contributions from s 
and u but have now the problem of the region ( > x > described in connection with the 
form factor. What we will do in this paper is to use LCWFs to calculate the contribution of 



5 Note however that in those cases there are factorisation theorems stating that the handbag diagrams 
arc dominant when the hard scale becomes infinitely large. In the present case we have a less strong 
situation of factorisation since for infinitely large —t the hard scattering mechanism of Q dominates over 
the soft overlap or handbag contribution. 

6 A more careful discussion is needed in the case where s ^> Q 2 , which we shall not consider here. 
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the lowest Fock state components to skewed parton distributions in the region 1 > x > (. 
We are thus not able to predict the amplitude of the DVCS process but can give a part of 
the nonperturbative input needed to calculate it, which furthermore is process independent 
and also occurs e.g. in exclusive meson production at large Q 2 and small t [Qj. 

It should be noted that even if we were able to express the full DVCS amplitude through 
the overlap of LCWFs we could not hope to evaluate the amplitude from the lowest Fock 
states only. In the case of the elastic form factor, where we do have an overlap formula, we 
know that all Fock states become important as one goes to low —t and it seems reasonable 
to expect the same for DVCS, where —t is always small by definition. Similarly the usual 
parton distributions, where we have an overlap formula in the full range < x < 1, can 
be well described by the first few Fock states down to some finite value of x, but at some 
point higher Fock states will become essential. The same holds a fortiori for skewed parton 
distributions as we shall see in Sect. [| 

3.2 Cat's ears diagrams in Compton scattering 

So far we have only considered soft overlap contributions with only one active parton, 
which is subsequently hit by the two photons. As already remarked they have the topol- 
ogy of handbag diagrams, i.e. they factorise into a parton-photon scattering and a soft 
subamplitude with two hadron and two parton lines, which we want to describe in terms 
of hadron LCWFs. There are other overlap contributions with two active partons, each 
coupling to one photon; they have the topology of so called cat's ears diagrams. One can 
see that in the large angle region as well as for DVCS one cannot avoid large virtualities or 
intrinsic transverse momenta occurring somewhere in these diagrams, so that we no longer 
deal with a soft overlap. Working with soft hadron wave functions one must then add at 
least one hard gluon in the diagrams. 

That in DVCS cat's ears diagrams become unimportant in the large-Q 2 limit is part of 
the factorisation theorem for that process. In the large angle region it is interesting to note 
that the diagrams where there is just one hard gluon exchanged between the two active 
partons consist of a hard scattering subprocess involving two parton lines (corresponding to 
the diagrams for Compton scattering off a meson in the hard scattering mechanism of 0]), 
and a number of spectator partons which as in the soft overlap (handbag) diagrams must 
be wee partons. It is reasonable to assume that such "hybrid" diagrams give contributions 
to the amplitude whose order of magnitude is between the pure soft overlap and the pure 
hard scattering contributions: compared with the latter they have less hard gluons (and 
thus hard propagators and coupling constants), but in contrast to the pure soft overlap 
diagrams they require N — 2 instead of N — 1 wee partons in an iV-particle Fock state, 
which is less restrictive for the hadron wave functions. 

3.3 Soft overlap contributions to other processes 

Having discussed in detail the conditions necessary to express soft overlap contributions in 
terms of LCWFs for spacelike elastic form factors and for Compton scattering we wish to 



13 



make some remarks on other processes: 



3.3.1 Meson production Y'p — > Mp 

Let us first see what happens if in the overlap diagrams for Compton scattering we replace 
the outgoing photon with a meson M = p, <fi, ir, K, . . . , and the pointlike photon-quark 
coupling with the qq Bethe-Salpeter wave function of M. In the discussion of our toy model 
we have seen that the loop integration over k~ gives a sum over residues, where each term 
corresponds to a simple pole in the /c~-plane and can be written as the product of two 
LCWFs (of the two external particles that "share" the parton which is on its mass pole). 
This is not the structure we would need for an expression in terms of three LCWFs, two 
for the incident and scattered proton and one for the meson. If and how such a structure 
can be obtained requires further investigation which goes beyond the scope of this paper. 

From our discussion of Compton scattering it is however clear that in the region of 
large s, —t, —u there is no soft overlap, because if the partons in the protons are all to be 
soft then there is a parton with large virtuality s or u, which now couples to the meson. 

In the region of small —t but large Q 2 and s the situation is different. First we remember 
that in this case it has been shown |24]] that for longitudinal photon polarisation and in the 
large-Q 2 limit the process factorises into a soft amplitude involving the two protons and 
two partons, the soft transition from a gg-pair to the meson, and a hard photon-parton 
scattering with at least one hard gluon exchange. A soft overlap contribution competing 
with this mechanism is possible when the quark line that directly goes from the meson to 
the soft proton amplitude is a wee parton: then one can take out the gluon from the hard 
scattering diagrams without any parton line going far off shell .[] As mentioned above it is 
not clear whether such contributions can be written in terms of LCWFs for the protons 
and the meson. Likewise it remains to be investigated whether it can be expressed in terms 
of the LCWF of the meson and a skewed parton distribution in the proton, the latter being 
obtained from the parton-proton amplitude by an integration over k~ in a similar way as 



LCWFs are obtained from Bethe-Salpeter wave functions p5| 



3.3.2 Timelike processes 

Crossing relates the timelike (7* — > pp) to the spacelike form factor (^y*p — > p), and the 
production of pp in a two-photon collision to Compton scattering; the diagrams for the 
timelike processes are obtained from those in Fig. [I] by a rotation of 90° counterclockwise. 
Using our toy model one easily sees that like their spacelike counterparts these processes 
admit soft overlap contributions. They can however not be expressed in terms of LCWFs: 
the parton line shared by the proton and antiproton cannot correspond to an incoming 
parton for both p and p as it would have to be in LCWFs, except for the the point where 
its plus momentum is strictly zero. This holds in any reference frame so that knowledge of 

7 Such end point configurations are indeed the reason why factorisation cannot be established in the 
case of transverse photon polarisation. 
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the LCWFs is not sufficient to evaluate the soft overlap contributions to these processes.^ 



4 Large angle Compton scattering with the handbag 



4.1 Calculation of the handbag diagrams 

The calculation of the handbag diagrams for real or virtual Compton scattering at large s, 
—t and —u can be done using the methods that are well known for usual DIS and for DVCS. 
At some points it presents however additional complications which we shall now discuss. 
For simplicity we work in the frame of Sect. |2.2| where ( = 0, although our derivation 
can be done in other frames as well. Our starting point is the expression of the Compton 
amplitude in terms of a soft proton matrix element and a hard parton-photon scattering: 



+ (p'\ T^j aa (z) W(0) \p) H a/3 (-k, —k') 



where 



H a0 (k', k) 



( fc + <?) • 7 

e -7- 



(k + q) 2 + ie 



e ■ 7 + e ■ 7 



{k'-q)-l 
{k' — q) 2 + ie 



e -7 



(23) 



(24) 



a/3 



is the tree level expression for the hard scattering, with polarisation vectors e and e' for 
the incoming and outgoing photon. The sum is over quark flavours a, e a being the electric 
charge of quark a in units of the positron charge e. The first term in (^3j) corresponds to 
the case where the incoming parton k in the hard subprocess is a quark, the second term 
corresponds to an incoming antiquark. For ease of writing we do not display the spin labels 
for the proton states here and in the following. 

Using that the photon-parton scattering is dominated by a large scale we now neglect 
the variation of the transverse and minus components of k and k! in H, where we replace 
them with momentum vectors that are on shell and lie in the scattering plane, namely with 



t/A 



■2k + " "2 Aj 



k' 



t/A 



2A-+ ' 2 Aj 



(25) 



respectively. The integrations over k~ and in (^) can then be performed explicitly, 
leaving us with an integral J dk + J dz~ and forcing the relative distance of fields in the 
matrix elements on the light cone, z — [0, z~, OjJ. After this the time ordering of the fields 
can also be dropped |25 . 



At this point one might be tempted to proceed as in standard DIS (or in DVCS) and 
decompose H on the Dirac matrices r y p and 7 p 7s- This leads to the Fourier transforms of 



8 Again it might be possible to find an expression of the overlap along the lines mentioned in our 
footnote but this will not be pursued here. 
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the nonlocal matrix elements {p'\ ip a (0) 7 P ijj a {z) \p), (p'\ ^a(O) 7 p 7s ip a (z) \p), and the corre- 
sponding ones with the arguments and z interchanged. In DIS or DVCS, where only the 
plus components of the proton momenta are large, one has that only the plus components 
of the currents give a leading contribution in the limit of large Q 2 . Now however we have 
a large scattering angle, and the proton momenta have large plus, minus and transverse 
components, so that it does not follow from kinematic considerations that the plus compo- 
nent of, say, / dz~ e lk+z (p'\ V'a(O) 7 P ip a ( z ) \p) is large compared to its minus or transverse 
components and thus dominates in the Compton amplitude. 

To show that the plus components indeed dominate also in large angle scattering we use 
that the proton-parton amplitudes described by the soft matrix elements can be written as 
the amplitude for a proton with momentum p emitting the active parton with momentum 
k and a number of on-shell spectators times the corresponding conjugated amplitude for 
momenta p' and k', summed over all spectator configurations; this just corresponds to 
inserting a complete set of intermediate states between the quark and antiquark fields in 
the matrix elements. We note that for small k 2 , k' 2 and small intrinsic transverse parton 
momenta, V 2 1 _ i jx i ^ A 2 , one cannot form large kinematical invariants at the hadron-parton 
vertices.^ 

For each of the proton-parton vertices we now go to a frame where the momentum k or 
k' has a zero transverse (and thus also a zero minus) component, performing a transverse 
boost as described in Sect. |2T4] . Considering for definiteness the case where the parton 
coming out of the proton is a quark we write in this frame 

i[)(z) = ^7~7 + ^0) + ^7 + 7~ ip( z ) 

= E [«(*. A ) ^7 + ^)) + iMk A) (u(k, XMz))} (26) 

A 

with a sum over helicities A/2 = ±1/2. We can now argue that in the matrix element 
of (|26|) between the incoming proton and the spectator system the term with u / ~f + ifj(z) 
dominates over the one with uip(z) because at the vertex we have a large plus component 
but no large invariant, and thus retain only the first term in the decomposition (p6|) .f^] 
Now we use that the boost ([13|) to the frame where k has vanishing transverse and minus 
components leaves the plus component of any vector unchanged so that fl2"6|) also holds in 
the overall symmetric frame of Sect. |2]^. Repeating our argument for the antiquark field 
we arrive at the replacement 

/ 1 ^ 2 



MO) M*) - Up) £ (mi + u{k\ A')) («(*, A) 7 + ^(^)) x u a (k', A') up(k, A) (27) 



A, A' 



and an analogous one involving antiquark spinors for ip a (z) ip a i (0). In (|23D the hard scatter- 
ing kernels are then multiplied with the spinors for on-shell (anti) quarks, which guarantees 

9 The situation is special for small momentum fraction x of the active parton, when Fock states with 
large N are important; a case we do not consider here. 

10 We note that this corresponds to the "good" component of the Dirac field in the context of light cone 
quantisation p3. 
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electromagnetic gauge invariance of our result. Note that the full expression (^3]) need not 
be gauge invariant since the handbag diagrams are not the complete set of diagrams for 
our process. 

To further simplify the hadronic matrix elements we use that the hard scattering, where 
of course we neglect quark masses, conserves the parton helicity: A' = A. In a suitable 
convention for massless spinors one has u(k, A) = —v(k, —A) and arg X)^ + u{k 1 A)] = 1 
for any on-shell momenta k, k', so that we can multiply (p7|) with 



w(fc / ,A)7+M(fc,A) 
1 = 2^ • (28) 

With u(k, A) u(k, A) = k • 7 (1 — A7 5 )/2 and analogous expressions for k' and for antiquark 
spinors we obtain after a little algebra 



X 



(p'\ MO) 7 + ^^ \P) X)H{%, k)u(k, A) 



+ (p'l M*~) 7 + ^^ ^a(O) \p) v{k, X)H{-k, -k')v{k', A) 



(29) 



where we write ip(z~) as a shorthand notation for i/j(z) with z — [0,z~ , 0j_]. We thus find 
that the plus component of the nonlocal currents dominates as we have anticipated in our 
footnote |3], and that the operators in the matrix elements are in fact the same as those of 
the leading- twist parton distributions occurring in DIS or DVCS. 

We now must discuss what to take for k + in the hard scattering. As shown in Sect. p. 1.2 



the requirement to have no hard partons directly coupling to the protons forces the active 
partons k and kl to have small intrinsic transverse momenta in their parent hadrons and a 
momentum fraction x close to one when —t is large. This corresponds to the approximation 
(p5|) with k + = p + we will make in the hard scattering factors, i.e. the expressions after 
the proton matrix elements in (|29|). Some degree of arbitrariness is associated with the 



global factor l/(2/c + ) in (p9|), which has its origin in (|2q) , (|27D and for which we choose 
to keep k + = xp + . Admittedly there is no clear-cut way to associate it to either the hard 
scattering, where we set x = 1, or the soft matrix elements, where setting x = 1 would not 
even make sense since for x strictly at its end point our proton LCWFs are zero. 

Making use of the charge conjugation properties of Dirac matrices and spinors in or- 
der to rewrite the term corresponding to antiquark-photon scattering we obtain our final 
expression for the handbag diagrams in large angle Compton scattering, 



A = - A Y.u{k\\)H{k\k)u{k,\)Y,M 2 I — f^e lxp+z ' 
4 *y a Jo x J 2it 

x [(p'\ Mo) 7 + M*~) - Mz~) 7 + Mo) \p) 

+ x(p'\M0)i + i 5 MM + MMi + i 5 Mo)\p)] , (30) 
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with (E3J) and with ( 25 ) for k + = p + . We note that the Fourier transformed matrix elements 
in ( p0[) are skewed parton distributions at ( = and large —t, as was already remarked 
in Jl^|. In (3C) we have incorporated their support property x < 1, cf. 0, 25 1. Following 



Radyushkin |L7j we introduce a form factor decomposition 
" l dx , f dz~ 



„ Jo x J 



2tt 



3 i xp^z 



n «/ U *JU J 



dz 
~2^ 



3 « xp^z 



(p'| Va(0) 7 + - 7 + ^a(O) \ P ) 

= R v {t) u{p') 7 + n(p) + R T (t) ±- u{p') a + »A u u{p) , 

2m 

(P'\i'a(0)l + 75lpa(z~) +4'a(z~)'J + 'J 5 1p a (0) \p) 

= R A (t)u(p') 1 + l5 u(p) +R P (t)^-u(p') l5 u(p) (31) 

2m 



for the x-integrals over these skewed distributions. Ry, Rt, Ra and Rp are new form 
factors specific to Compton scattering; note that Rp does not contribute to the Compton 
amplitude in our symmetric frame with A + = 0. 

One may ask how to improve on the approximation (|25| ) with k + = p + when calculating 
the hard scattering. There will be corrections due to the facts that in the hard scattering 

1. x is not strictly one, 

2. the intrinsic transverse momenta of the partons k, k' are nonzero, and 

3. the virtualities k 2 , k' 2 are not zero. 

The order of magnitude of all these corrections is controlled by the parameter A 2 / (— t) as 



discussed in Sect. |3.1.2| . Note that in order to express the amplitude in terms of LCWFs or 
of the light cone matrix elements in (^) it was essential to neglect the ^"-dependence of the 
hard scattering. The inclusion of off-shell corrections (point 1) would thus necessitate an 
extension of the framework we are using here. We emphasise that the on-shell condition in 
the hard scattering is our guarantee to obtain a gauge invariant result; "exactly" evaluating 
the handbag diagrams would only have a limited sense since a part of the corrections to 
(f25|) will break gauge invariance and be cancelled by other diagrams. Furthermore points 1, 
2 and 3 are kinematically related: from k 2 — k' 2 = 2A^ • (kj_ + Aj_/2) in our symmetric 
frame we see that if we insist on taking on-shell partons in the hard scattering then we 
must fix kj_ = — Aj_/2 (as we did in (p5[)), which forbids us to evaluate the effect from 
the variation of kj_ in the hard scattering kernel. We also see that for x ^ 1 the choice 
h± = —A.±/2 no longer corresponds to zero intrinsic transverse momenta k± + xA.±/2 and 
kj_ + (2 — x) Aj_/2 of k and k' in their parent hadrons. 

Compared with fixing k^ the approximation x — 1 in the hard scattering presents the 
particularity that x is taken at its kinematical end point; the soft part of the process can 
only select x around some value smaller than 1. Moreover, we find that with our ansatz for 
the LCWFs (Sect. |5]) both the ^-integrals in (|3lD and the corresponding one for the elastic 
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form factor are dominated by values of x not very close to 1 for —t between, say, 5 and 20 
GeV 2 , with the peaks of the integrands being of order 0.45 to 0.75. The reason is that with 
our wave functions the end point x — 1 is rather strongly suppressed in the integrands of 
fl68|) , ([T3|) by a third power (1 — x) 3 , cf. (^|), (0), (^1), and that the suppression of large 
k. 2 ± /x in the LCWFs is only effective for values clearly larger than 1 GeV 2 . It turns out 
that the factor 1 jx in (|31|) does not significantly shift the values of x where the integrand 
has its maximum, but rather increases the height of the peak. 

One might think of only dropping the approximation x = 1 then, but allowing x to 
be different from 1 in the hard scattering would lead to serious problems: in the case 
of a real incident photon for instance one easily calculates that for = — Aj_/2 and 
x = r] = (y/s — y/—u)/ (y/s + \J — u) one has u = 0. It would however be mistaken to treat 
this as a pole in the hard scattering ([24]) which gives an imaginary part to the scattering 
amplitude. We must remember from our discussion of the ^"-integration in Sect. |3.1.2| 
that we have already neglected certain terms where u has a pole. Retaining others by 
allowing x to range from to 1 in the hard scattering is then inconsistent and would 
give misleading results. What happens in this example is that the factorisation into a 
hard scattering and a soft proton matrix element breaks down for x not sufficiently large. 
Keeping x — 1 fixed in u(k')H(k', k)u(k) is thus related to our approximation of factorising 
the soft overlap contribution to Compton scattering into a hard parton-photon scattering 
and a soft proton matrix element. 

The fact that in our numerical applications the hadron wave functions are probed at 
intermediate rather than very large x means on one hand that our results are not too 
sensitive to the precise behaviour of the LCWFs near x — 1, and also not to a possible 
Sudakov suppression (cf. [16j for comments on these points in the case of the elastic form 
factor). On the other hand our approximation x = 1 in the hard scattering of the Compton 
process has only a limited accuracy for —t not very large. 

We finally also neglect the proton mass when relating s and u to the external variables. 
Comparing (|]) with (|25| ) at x = 1 we see that this means k w p and k' ~ p' so that we 
have s ~ s and u ~ u. Corrections to this will be of relative order m 2 /(—t) and thus of 
the same size as other terms we do not control. 



4.2 Proton spin 

We have already remarked that the hard scattering subprocess does not change the helicity 
of the active parton (the same holds for the quark-photon coupling in the elastic form 
factor). As the helicities of the spectators do not change either a change in the proton 
helicity implies that for at least one of the incident or scattered proton the parton helicities 
do not add up to the hadron helicity. In other words the calculation of proton spin flip 
amplitudes requires to take into account LCWFs with nonzero orbital angular momentum 
L3 of the partons in a detailed manner^ this will not be attempted in the present work. 
For the lowest, three quark Fock state we only take a wave function with zero L3, which has 

11 In this respect one has the same situation as in the hard scattering formalism p6fl. 
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been constrained by several physical observables in and do not endeavour to model 
wave functions with L 3 7^ 0. For higher Fock states, which for sufficiently large —t provide 
only a correction to the three-quark contribution in Compton scattering and the elastic 
form factor, we will not specify how the orbital angular momenta between the various 
partons are explicitly coupled; describing such detailed effects is not within the scope of 
this paper. 

Due to its finite mass the helicity of a proton depends of course on the choice of reference 
frame. Taking the incident proton for definiteness we can express this dependence in a 
covariant way using its spin four-vector s. In the hadron frame of Sect. |2.4j , where p 
has zero transverse momentum, the spin vector for a state of definite helicity is a linear 
combination of p and the vector v' — [0, 1,0 J, which is unchanged by the boost to the 
overall symmetric frame. This choice of spin quantisation axis is natural in our context 
of LCWFs, which are defined with respect to the same vector v' through the integration 
over the minus components of parton momenta. A corresponding argument holds for the 
scattered proton, with the same vector v'. We find that in our symmetric frame with ( = 
the helicity flip amplitudes are only due to Rt, which we will therefore not be able to 
model here, whereas the helicity conserving ones go with Ry and Ra- The same holds for 
the elastic form factors: F 2 changes helicity and Fi does not, and we will only calculate 
Fx. 

We know from experiment that in the transition 7*p — > p proton helicity flip becomes 
small compared with no flip for large enough —t, so that neglecting the former can be 
justified as an approximation. The measured difference between the Dirac form factor 
Fx and the magnetic Sachs form factor Gm = F 1 + F 2 at a given —t shows the degree 
of accuracy of neglecting spin flip contributions, and it is reasonable to assume that the 
situation will be similar for the new form factors fl3"T|). 

4.3 The hard scattering 

We now give the hard scattering amplitudes 

H x> ^ = u(k',\)H(k',k)u(k\) (32) 

where \x and // respectively denote the helicity of the initial and final state photon. For 
virtual Compton scattering the initial photon helicity depends on the reference frame and 
we choose to define it in the photon-proton cm., i.e with respect to the p-q axis: our 
symmetric ( = frame is adapted to discuss the physics of our reaction mechanism, but 
7*-polarisations defined in the cm. are well suited for the consideration of azimuthal asym- 
metries we shall briefly mention below, apart from being a standard choice that facilitates 
comparison with other work. With our approximation k ~ p, k' ~ p 1 the photon-proton 
cm. is identical to the cm. of the hard subprocess q(k) 7* (9) — > q(k') 7(9') • m our phase 
convention, where axg[u(k', X) , ~f + u(k, A)] = 1, we explicitly find 

nr s + q 2 

— 2 4/ , H 
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S + Q 2 y/—SU ' 

W+,o- = -2 — -7^ y/—2t , ^+,0+ = 0, (33) 
s + Q 2 

with the kernels for A = — 1 given by parity invariance as Hx,^ = (— l) M-/i W_a, -h-^'- 

With ( p0| ) and ( |33"D we have all necessary ingredients to calculate the cross section in 
terms of the form factors Ry, Ra (and Rt which we will neglect in this work). We present 
a numerical study of real Compton scattering in Sect. |8|. Virtual Compton scattering 
is measured in electroproduction, ep —>■ epj, where it interferes with the Bethe-Heitler 
process, i.e. the emission of the final state photon from the lepton, and its detailed study 
shall not be attempted here. The results in the handbag mechanism have however some 
general features, both for real and virtual initial photons, which we discuss now. 

The first point is that the photon-proton amplitude comes out as purely real: the form 
factors Ry, Ra, Rt are real due to time reversal invariance, and the hard scattering kernel 
does not have an imaginary part because the corresponding diagrams cannot be cut with 
s and u being far off-shell; such cuts only arise at the level of a s -corrections to the photon- 
parton scattering. In the hard scattering mechanism [[| the situation is very different: 
there one has cuts already to leading order in a s , which lead to nontrivial phases in the 
scattering amplitude. This may offer a valuable tool to distinguish experimentally which 
reaction mechanism is at work: in ep —>■ epj with longitudinally polarised lepton beams 
the beam polarisation asymmetry is proportional to the imaginary parts of the 7*p — > 7p 
helicity amplitudes, with the Bethe-Heitler amplitude being purely real. In the handbag 
mechanism this polarisation asymmetry is then predicted to be small, arising only at the 
level of loop corrections, while in the hard scattering mechanism it can have a substantial 
value. This was for instance shown in [[^J, where virtual Compton scattering was studied 
within the hard scattering approximation using a quark-diquark wave function for the 
proton. 

A second remarkable feature of ( |3"3"D is the dependence on the photon helicities: transi- 
tions between positive and negative helicities are forbidden for real photons and suppressed 
by Q 2 /(s + Q 2 ) if the photon virtuality is small compared with s.Q This helicity selection 
rule could be tested in real Compton scattering with linearly polarised incident photons: 
it leads to the absence of a dependence of the cross section on the azimuth $ between the 
plane of photon polarisation and the scattering plane; nonzero photon helicity flip ampli- 
tudes will in general give a cos2$-contribution to the differential cross section. For finite 
but not very large Q 2 the situation is more complicated, because the helicity flip amplitude 
is only suppressed and not zero, and because the process interferes with Bethe-Heitler. 



Whether this still holds at the level of a s -corrections would need further investigation. 
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5 The Fock state wave functions 



The valence Fock state of the nucleon has been investigated in some detail in Ref. [16 . 



The explicit form of the corresponding wave function has been extracted from a fit to the 



valence quark distribution functions derived in Ref. |28[ an d to the Dirac form factor of 
the proton assuming dominance of the soft overlap contribution. This is just the physics 
we are interested in here; therefore we take over the results of Ref. |16| as a starting point. 
The wave function proposed in Ref. |16|] has also been shown to work successfully for J/ip 
decays into proton-antiproton pairs, a process that is well under control of perturbative 
physics in contrast to, for instance, the form factors in the experimentally accessible region 
of momentum transfer. In the subsequent sections we will test that wave function in further 
observables, namely in Compton scattering and in the polarised parton distributions. We 
will even go a step further than in Ref. [16| and explore the next two higher Fock states 



consisting of four and five partons in order to determine their gross features. Moreover, 
we are going to investigate the global effect of all Fock states in an approximate way. As 
has been shown recently by Radyushkin |17|, one can then directly relate the parton dis- 



tributions controlling deeply inelastic lepton-nucleon scattering with exclusive observables 
such as form factors or the Compton cross section, without assuming an explicit form of 
the distribution amplitudes. 

For the reader's convenience we will start with a brief description of the properties of 
the LCWF for the proton's valence Fock state derived in |16|. According to Sotiropoulos 



and Sterman [2~S| the valence Fock state of a proton with momentum p and positive helicity 



can be written as 

|P,p,+;qqq) = J [dx] 3 [d 2 k ± ] 3 {^ 123 M+-+ + ^ 213 M-++ 

-(#132 + ^23l)M++-} , (34) 

with plane wave exponentials and the colour wave functions omitted here and in the fol- 
lowing. The integration measures in Eq. ([34]) are defined by 

N i N 

[dx] N = J] dx, 6(1 - 5><) , [d 2 k J,v = 3 I] 5(2) (E k ^*) • (35) 

i=l i l iC)7r ) i=l i 

The quark i is characterised by its plus momentum kf = Xip + , its transverse momentum 
kjj with respect to the proton's momentum, and by its helicity Aj. A three-quark state is 
then given by 

M Xl x 2 x 3 = ; I u; xip + , k ±1 , A x ) | u; x 2 p + , k ±2 , A 2 ) | d; x 3 p + , k_L 3 , A 3 ) (36) 

y/XiX 2 X 3 

with a normalisation 

(q; k' ±i , A- 1 q; Xip + , k ±i , X t ) = 2xip + (27r) 3 5 A / A .<5 (x[p + - x t p + ) 5(k' ±i - k ±i ) . (37) 
A neutron state is obtained by the exchange u <-> d. 
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We only consider the part of the wave function with zero orbital angular momentum L 3 
along the 3-axis, so that the quark helicities sum up to the proton's helicity. As has been 
demonstrated in Ref. |3(J Eq. fll^f ) is the most general ansatz for the L3 = projection 
of the three-quark proton wave function: From the permutation symmetry between the 
two u-quarks and from the requirement that the three quarks have to be coupled in an 
isospin 1/2 state it follows that there is only one independent scalar wave function, which 
for convenience is parametrised as 

/ 3 

123 fat, k_u) = *(ari, x 2 , x 3 ; k ±1 , k ±2 , k ±3 ) = — = <pi 23 (xi) ^ 3 (x u k_u) (38) 

with the normalisation conditions 

/ [dx} 3 <f> m { Xi ) = 1 , J [d 2 k ± ]jv Q N (xi, k ±i ) = 1 . (39) 

/3 plays the role of the nucleon wave function at the origin of coordinate space and 
4>i23(xi) = 4>{x\,X2,x 3 ) is the nucleon's valence distribution amplitude. Both quantities 
depend on a factorisation scale fip and are subject to evolution. Expanding 4>i 23 (xi) as 

(f>123(Xi, Hf) = 4>As(Xi) 1 + B 1 (fM F ) 0i 2 3 (Xi) + B 2 (fi F ) 0123 ( X i) + •• • , ( 40 ) 

where 4>As(xi) = I2OX1X2X3 is the asymptotic distribution amplitude and 4>\ 2 3(xi), <fii 23 {xi), 
etc. are the eigenfunctions of the evolution kernel H, one has 



2/(3/3o) / w /A \ \ 7n//3o 

(41) 



with p = 11 — 2n//3, 71 = 20/9, 72 = 8/3, etc. In Ref. |L6| it has been shown that it 



is sufficient to retain only the first two terms in the expansion ( j4D[) . They are taken as 
Bi(no) = 3/4 and B 2 (fio) = 1/4 at a factorisation scale of /xo = 1 GeV. At this scale one 
then has the simple form 

(pn 3 (xi) = 60 X1X2X3 (1 + 3xi) (42) 

for the valence distribution amplitude. 

When calculating the overlap contributions to the elastic form factor and to large angle 
Compton scattering in Sect. |7| and |8| we will use the distribution amplitude at a factorisation 
scale n F = —t given by the momentum transfer to the nucleon.^] The parton distributions 
in Sect. ^ and || will be calculated and compared with the parametrisations from global 
fits at our starting scale /j,q = 1 GeV 2 . 

The transverse momentum dependence of the wave function is contained in the function 
f2/v- A simple symmetric Gaussian parametrisation, 



n t ^ \ 16 ^ a NJ 
U N {Xi,k ± i) = exp 

XiX 2 . . .x N 



N 12 

a — 1 %i 



(43) 



13 



For the higher Fock state LCWFs to be discussed below the evolution will be neglected. 
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suffices and meets various theoretical requirements, see for instance pi [31], |32| and our 



remark following Eq. ([22]). This ansatz keeps the model simple and allows one to carry 
through the kj_-integrations analytically. Note that the ansatz fl38|), ( f43|) represents a soft 
wave function, i.e. the full wave function where the perturbative tail with its power-like 
decrease is removed [|J. Integrating in Eq. (^) to infinity instead of to a cut-off scale 
given by the hard scale in a process introduces only a small negligible error. 

The values of the normalisation f 3 and the transverse size parameter a 3 have been 
determined in [[Toll as 



f 3 = 6.64 ■ 1(T 3 GeV 2 , a 3 = 0.75 GeV" 1 (44) 

at the scale of reference \ip = Ho- With these parameters the valence Fock state wave 
function has a value of 0.17 for its probability and a value of 411 MeV for the rms transverse 
momentum. The valence Fock state thus appears to be rather compact, with a radius of 
only about a half of the charge radius. For further discussion of the properties of the 
valence Fock state wave function see []TB] . 



With the valence Fock state fully specified we can now turn to the higher ones. Ex- 
plicitly we only consider the Fock states with an additional gluon (N = 4) and with an 
additional sea quark-antiquark pair (N = 5). Due to parity conservation both require one 
unit of orbital angular momentum. One therefore encounters many different possibilities 
of coupling the various partons in a nucleon, each coming with a new wave function. It 
seems plausible to assume that the effect of the orbital angular momentum is averaged out 
in the sum over all different coupling possibilities.^ 

With this proviso in mind we take 



P,p,+;qqqg) = y"[dx] 4 [d 2 kj 4 ^1234(2^, ky 1 



x [M++- - M+-+] — L I g; x 4 p + , ku, A 4 ) (45) 



as a representative of all N = 4 Fock states, with the gluon state | g; x 4 p + , kj_ 4 , A4) nor- 
malised as in ( |3"TP and 

*1234(Zi, k_u) = — ^= 01234(^) ^4(^1, k±i) . (46) 

For the distribution amplitude of this Fock state we take (at the scale /zo) 

9! 

4>1234{Xi) = — XiX 2 X 3 X 4 (1 +3xi) , (47) 



14 In principle there is no difficulty in treating all possibilities explicitly. Each of them is described by 
an appropriate covariant spin wave function |3^, |34|] that is proportional to K ■ 7, where is the relative 
momentum of two clusters of partons. These K ■ 7-terms, representing the orbital angular momentum 
between the two clusters, give rise to an additional factor ~ K'^ in the expressions for observables like 
the overlap integral for the nucleon form factor. 
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i.e. the distribution amplitude has the asymptotic form multiplied by an asymmetry factor 
of the same type as in the distribution amplitude for N = 3. The spin-isospin coupling of 
the valence quarks requires a distribution amplitude that is symmetric under the exchange 
2^3. The gluon is supposed to couple with the orbital angular momentum in a spin zero 
state. Thus the ansatz ( |45| ) satisfies the minimal requirement that the partons of this Fock 
state are coupled in a spin-isospin 1/2 state. 

For the N = 5 Fock state we assume a sea that is colourless, SU(3) flavour sym- 
metric and coupled to total spin zero.0 The generalisation to a more complicated sea 
is straightforward, requiring flavour-dependent wave functions which may also have ad- 
ditional asymmetries in their Xi -dependence, but in order to keep the model simple we 
refrain from introducing such wave functions. With our simple ansatz the valence quarks 
are in a totally symmetric state in flavour-spin-momentum-space, just as the valence Fock 
state itself, and the valence sector of the N = 5 Fock state therefore exhibits the same 
structure as fl54"|). Assuming its wave function to equal that of the valence Fock state we 
make the ansatz 

| P,p,+;qqqqq) = y[dx] 5 [d 2 k ± ] 5 

X {^12345 + ^21345 M- ++ ~ (^13245 + *2314 5 ) M ++_ } 

x {|q; x 4 p + ,k ±4l +) |q; x 5 p + ,k ±5} -} 

q=u,d,s V 4 5 

- I q; x 4 p + , k_L4, -) | q; x 5 p + , k_L 5 , +) } (48) 

with 

fs 

^12345(2^, kj_i) = — 012345 fai) ^5(2^, k ±i ) (49) 
for the wave function and 

10! 

012345 (#i) = — X 1 X 2 X 3 X i X 5 (1 + 3x x ) (50) 

Id 

for the distribution amplitude at scale /io- The symmetrisation between the sea and valence 
quarks required by the Pauli principle is ignored here. We argue that it cannot play a major 
role because of the fairly large spatial separations between sea and valence quarks: the sea 
quarks have to build up the full charge radius of the nucleon while the valence quarks form 
a compact core. 

Admittedly our parametrisations of the higher Fock states are oversimplified. For the 
physical processes and in the kinematical region of interest here they give however only 
small contributions compared with the valence Fock state, and to investigate these correc- 
tions we deem our ansatz to be sufficiently accurate. 

15 The N = 5 Fock state with two gluons in it is discarded since its contribution to physical quantities is 
highly suppressed in the kinematical region of interest, cf. our remark after (p^) below. If in the following 
we talk about higher Fock states (N > 5), this particular Fock state is understood to be included. 
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We finally give an integral which will appear in our overlap formulae, namely 



I N (xi, C, A* ) = J [d 2 k±] N Q N (x'i, k' ±i ) Q N (xi, kj 



(51) 



where the relation between the primed and unprimed variables is given by (16) and for £ = 
also by flI5|). An appropriate tilde, hat or breve upon the variables Xj,k^j is understood 
depending on the case. For our Gaussian k^-dependence the integral (|51|) evaluates to 



Pn 



x i 



X' 



-T^(x„C;Ai) 



(52) 



N 



where 



N-2 



r N (x h (;Al 



2{ Xj - C) 



CJ fa - Q + Xj (i - c)'< 



exp 



~ a% Aj (1 - Xj ) 
( Xj -C) + x 3 (l-0 2 \ 



and 



Pn 



2„2 sJV-1 



57T a 



A J i 



(53) 



(54) 



Notice that Jjy is a function of all momentum fractions Xj whereas T^v only depends on 
the fraction Xj of the active quark. 

Turning to a more generic notation we have that each Fock state is described by a sum 
of terms, each with its own momentum space wave function where f3 labels different 
spin-flavour combinations of the partons. On the basis of this notation the Fock state prob- 
abilities are given by P N = /[d^]jv[d 2 kj.]jv \^Np(%i, k_u 
we have 

435 27972 685125 

•T3 — ~Z~Z~7 p3 J 3 , -T4 — /?4 J4 ) ^5 — oc - P5 J| 



224^^ ' 275 ' " 4 352 

and with (|44[) we obtain P3 = 0.17 as already mentioned above. 



2 . For our parametrisations 

(55) 



'5 j 



6 Parton distributions 

As shown by Brodsky and Lepage H the contribution of an iV-particle Fock state to the 
distribution function for a parton of type a in the proton is generically given by 




where the sum j runs over all partons of type a. Summation over all Fock states leads to 
the full distribution functions 

q.(*) = E^O*)- (57) 

N 

Note that in our notation q a stands for the distributions of quarks, antiquarks or gluons. 
From the wave functions defined in Sect. |^ and with the help of ( |52"D for ( = and 
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Table 2: Coefficients for the Fock state contributions to the parton distribution functions 
according to Eq. (p8|). The powers m a (N) of (|61~D and (^) are also listed. 

Aj_ = one easily finds the individual contributions to the distribution functions from the 
N = 3, 4, 5 Fock states as a function of the parton momentum fraction x: 

qW( x ) = b[ N) P N x n " (1 - x) ma{N) [l + c<M (1 - x) + dW (1 - x) 2 ] , (58) 

where the coefficients c^' and are compiled in Tab. 0. 

As usual we define a valence quark distribution by q( N \x) = q^(x) — q^ N \x) for 
q = u, d. The sea is flavour symmetric in our simple model , hence 

u( 5 ) (x) = d (5) (x) = s (5) (x) = s (5) (z) . (59) 



With our particular ansatz (^8|) , (|50D we also have 

d {5 \x) = d^(x)/3. (60) 

One observes that all contributions appear in the form x n (1 — x) m times a polynomial 
in (1 — x) which is generated by the asymmetries in the distribution amplitudes, i.e. their 
departure from the asymptotic form. This holds for polynomial distribution amplitudes in 
general. The leading power m a (N) of (1 — x) in q^>(x) is generated by the corresponding 
asymptotic distribution amplitude; for quark distributions one has 



n, 



1, m q (N) = 2N + 2l g - 3 , (61) 



and for the gluon distribution 

n g = 3 , m g {N) = 2N + 2l g - 5 , (62) 

where l g is the number of gluons in the iV-particle Fock state. We see that higher Fock 
states generate higher powers m a (N). Summing over all Fock states the leading powers 
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of (1 — x) for valence quark, gluon and sea quark distributions come out as 3, 5 and 7, 
respectively. For the contributions from the N = 5 Fock state with three quarks and two 
gluons the leading powers are very high, m q = 11, m g = 9, which is why we do not consider 
it here. 

Our results for the valence parton distributions respect the usual counting rule be- 
haviour H. In other cases our results for the leading powers of 1 — x differ from those 
obtained from perturbative QCD arguments by Brodsky, Burkardt and Schmidt ^5[. This 
is not a contradiction since we are dealing with soft physics contributions. The perturbative 



results of Ref. |35[] manifest themselves only in a region 1 — e < x < 1 where the perturba- 
tive QCD contribution dominates over the soft contribution. To estimate e we remark that 
the overlap formulae fl56|) , (|67|) for parton distributions and elastic form factors are exact 
if one takes the full wave functions instead of their soft parts considered in this work ||. Q . 
Using the relations (|68|) or (]69[) between both types of quantities we obtain e ~ 1/ (— a\ i), 
where i is the momentum transfer in Fi p (t) at which the perturbative components of the 
wave functions start to dominate over the soft ones. For the wave function we consider 
here, —t is of the order of 500 GeV 2 JT5|. 



If for simplicity we take the transverse size parameters for the N = 3, 4, 5 Fock states 
to be equal, 

05 = 04 = 03, (63) 

then only one parameter remains free for each of the new Fock state wave functions, namely 
its probability (or the constants /at, cf. Eq. (^)). We fix these two parameters by fitting our 



gluon and antiquark distributions (|58"D to the Gliick-Reya-Vogt (GRV) parametrisation [|28 
at large x. A best fit is obtained for the values 

p 4 = p 5 = o.l , / 4 = 1.06 x 10" 4 GeV 3 , / 5 = 3.64 x 10~ 6 GeV 4 . (64) 

The results of the fit are shown in Fig. ^| and compared to the GRV parametrisation.^ The 
agreement with the distribution functions given in Ref. |36j is of similar quality at large 
x. All distribution functions of the proton are reproduced quite well down to x of about 
0.5, for the sea quark distribution even to lower values. We see how the first three Fock 
states control the large- x (x>0.5) behaviour of the distribution functions; certainly the 
situation could be improved by including even higher Fock states. We emphasise that the 
asymmetries in the distribution amplitudes play an important role: they push up u„ and 
diminish d„ at the same time, thus producing a ratio u v /d v of about five at large x while 
totally symmetric distribution amplitudes yield a ratio of only two, in sharp conflict with 
the GRV parametrisation. We also note that our ratio d v /u v tends to 1/14 in the limit 
x — > 1 and differs from the SU(6) result of 1/5 f37j. 

The spin dependent parton distributions allow another interesting test of our approach. 
These distributions measure the difference between the distributions of type-a partons with 
positive and negative helicity. In analogy to the unpolarised distribution discussed above 



16 At large x the 1998 GRV parametrisation is rather close to the 1995 version. We compare here with 
the LO parametrisation of the 1995 version. 
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parton distributions 

u v (x) 




parton distributions 

d v (x) 




parton distributions 

g(x) 




0.10 



parton distributions 

(ubar+dbar+sbar) 




Figure 3: Parton distributions obtained from the N = 3,4,5 Fock states (P3 = 0.17, 



P4 = P5 = 0.1). The model results are compared to the 1995 GRV LO parametrisation p8 



at a factorisation scale of 1 GeV. For the sea distributions we sum over the three flavours. 
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Table 3: Coefficients for the Fock state contributions to the spin-dependent parton distri- 
bution functions according to Eq. (|6q) . The powers m a (N) from (pi"|) are also given. 

we find within our model 

AqW(x) = A&jM P/v a; n » (1 - x) m ^ N) [l + Acf } (1 - x) + Adf > (1 - a;) 2 ] (65) 

with the coefficients listed in Tab. |3|. The powers n a and m a (N) are the same as the ones 
for unpolarised distributions, given by (|6~lf) . As a consequence of our simple assumptions 
that the gluons and sea quark pairs are unpolarised we have Ag^ 4 ^(x) = Aq^(x) = 0. Note 
also that Aq^'(x) oc q^'(x) at large x. While the constants of proportionality are close 
to unity for the valence u-quark distributions, they are negative or even zero for valence 
d-quarks. 

In Fig. |] we compare our predictions with the parametrisation proposed in Ref. J5S 



As we see, surprisingly good agreement is obtained in our simple model. There is also 
fair agreement with the polarised parton distributions determined in |311 at large x. The 



relative strength of Au v and Ad„ in that region reflects the spin structure of the valence 
Fock state and the asymmetry in its distribution amplitude. 



7 Form factors 



According to Drell and Yan [14[] the Dirac form factor can be represented as the overlap 
of LCWFs as 

Fi(t) = £ F{ N \t) (66) 

TV 

with individual Fock state contributions 

Fl N \t) = e a I [dxlivfd 2 ^]^^^,^)^^,^) , (67) 

a j 13 

where j runs over all partons of type a. We use our symmetric frame to evaluate the 
overlap, the primed and unprimed arguments in ( |5"T| ) are therefore related by (|T5p and we 
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Figure 4: Spin-dependent valence quark distributions Au v and Ad„. The model results 
are compared to the parametrisation of Ref. [SB . 



have = —t. Performing the k^-integrals for the N = 3, 4, 5 Fock states with the help 
of Eq. fl52]) , we arrive at 



Fl N n \t) 



j dx 


exp 


r 1 2 , 

2 a N l 


1 


— x~ 






j dx 


exp 


' l 2 + 
2~ a N l 


1 


— x~ 




X . 



{e u ui N \x) + e d di N \x)} , 
{e u di N \x) + e d ui N \x)} . 



(68) 



for the proton and neutron form factors. The appearance of the parton distributions 
here is a consequence of the fact that the integrand in their overlap representation (|5Bp is 
obtained from the one in (|67| ) by setting Ai = 0. Thus the k^-integrals only differ by the 
exponential factor of flSBI ) at ( = 0, which arises from the Gaussian k^-dependence of our 
wave functions. 

It is now suggestive to assume that the k^-dependence of all Fock state wave functions 
is given by the Gaussian (f£|) and to approximate all with a common transverse size 
parameter a. Summing over iV in (|66| ) then leads to a representation of form factors in 
terms of the valence quark distribution functions: 



F ln (t) 



j dx 


exp 


-a 2 t 
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1 — x' 


X 


j dx 


exp 


r 1 2 

-a 2 t 
2 


1 — x~ 


X . 



{e u M v (x) + e d d v (x)} , 
{e u d v (x) + e d M v (x)} , 



(69) 



a formula recently proposed by Radyushkin [I7|. Remarkably, inclusive observables are 
related to exclusive ones. The chief advantage of this formula is its independence from any 
explicit form of the distribution amplitudes. Of course a common value for the transverse 
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size parameter for all Fock states is unrealistic: as we saw before the valence Fock state 
is rather compact corresponding to about a half of the charge radius. Consequently the 
higher Fock states have to develop the full radius. For the purpose of evaluating the form 
factors from Eqs. ( |66]) and fl68[) we take (23 = 04 = as before and put as a simple ansatz 
a>N = 1.3 a 3 for N > 5, where the factor 1.3 is adjusted to the data for F\ p . A substantially 
larger factor would strongly suppress the higher Fock state contributions, a smaller one 
would lead to large contributions exceeding the form factor dataQ Then we set 

E qi N) (x) = q a (x)- E (70) 

N>5 JV=3,4,5 



where q a is taken from the GRV parametrisation |28[ and the three lowest Fock state 
contributions from our model. In this way we account for the sum of all Fock states in 
a phenomenological way. The results obtained in this manner are confronted to the data 
g in Fig. |. 




Figure 5: Electromagnetic form factor of the proton and neutron using the model parton 
distributions for the valence Fock state only, the N = 3,4,5 Fock states, and all Fock 
states on the basis of the GRV parametrisation at the factorisation scale 1 GeV ||28|| , cf. 
([7(f) . Data for Fx and Gm are taken from p0|, [41 



For large values of the momentum transfer our simple model agrees very well with 
the data, i.e. the dimensional counting behaviour is well mimicked by soft physics. Below 
about 10 GeV 2 the model is not perfect, deviations of the order of 20%, i.e. of the order 
of m 2 /(— t), are to be noticed. Such corrections are to be expected in our model, where 
proton spin-flip effects and orbital angular momentum in the wave functions are not taken 

17 We note at his point that in contrast to our ansatz a transverse size parameter a = 0.84 GeV -1 
common to all Fock states was used in [1171 . 
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into account as we discussed in see Sect. [4.2|. In reality spin flip effects are not very small 



as is indicated by the difference between the Dirac and magnetic form factors, F% and 
Gm-i see Fig. |5|. In view of these approximations we are satisfied with our results even 
in the range 5 GeV 2 < — t < 10 GeV 2 . We observe from Fig. [5] the dominance of the 
valence Fock state contributions. For —t > 10 GeV 2 all other Fock states contribute less 
than 20%; each individual Fock state provides only a small correction to the form factor. 
This can be regarded as a justification of the rough treatment of the N = 4 and 5 Fock 
states introduced in Sect. ^. We also remark that the parameters f 3 and a 3 of ([0]) used 
in this work have been obtained in [|16|] by requiring that the data for F± p be saturated 



by the soft overlap of the valence Fock state only. Given the uncertainties just discussed 
and our simplified treatment of the higher Fock states we think however that readjusting 
these parameters is not necessary here. As for the neutron form factor, we mentioned 
in Sect. || that totally symmetric wave functions lead to the relation u^ N \x) = 2d[ N \x), 
which according to Eq. ( p5[ ) would lead to a vanishing contribution to F ln . Hence the 
asymmetries in the LCWFs generate the neutron form factor. 

For wave functions of the type we are considering here the leading powers m q (N) of 
(1 — x) in the valence distributions o^\x) correspond to leading powers m q (N) + 1 of 
1/t in the asymptotic behaviour of the corresponding Fock state contribution to iq(t).[^] 
Hence for sufficiently large —t the valence Fock state dominates the form factor with only 
small corrections from the next Fock states. It is important to realise that this asymptotic 
behaviour of the overlap contributions does not set in before —t ~ 100 GeV 2 since the 
expansion of the integrals appearing in (]68f) into a power series in 1/t converges very 
slowly. We remark that the dominance of the soft overlap contribution is consistent with 
the strength of the perturbative contribution to the proton form factor, which drops as 
1/t 2 . As reported in Ref. [IB| the perturbative contribution evaluated from our valence Fock 



state wave function can be neglected for experimentally accessible momentum transfers. 
For —t larger than about 500 GeV 2 , however, the perturbative contribution will dominate 
since our overlap contribution asymptotically behaves as 1/t 4 . 

In analogy to the electromagnetic case we can also calculate the charged current axial 
form factor of the nucleon. The various contributions are now weighted by the quark 
helicities and isospin, leading to 



/ dx exp 


'1 2 1 — X~ 


o a Nt 




.2 x _ 



x {Au^(x) +2Au w (i) - Ad[ N \x) - 2 Ad {N \x)} , (71) 

Evaluating the axial form factor along the lines described for the electromagnetic case 
we find fair agreement with the dipole parametrisation of the admittedly low-t neutrino 
data Wl ■ 



18 The Drell-Yan result [jlij, for which the power m q of (1 — a;) in the valence quark distribution functions 
corresponds to a power (m q + l)/2 of 1/t in the form factor, is only obtained for wave functions factorising 
in x and k_L (i.e. for Q not depending on Xi). 
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8 Large angle Compton scattering 



Using our expressions (0), (|31~D for the handbag amplitude and neglecting the contribution 
from the proton spin flip form factor Rr, we obtain the cross section for real Compton 
scattering with unpolarised photons and protons as 



da 
dt 



2vra? 



u 
s 



- [Ry(t) + R A (t) 



us 



(R v (t) - R 2 A (t) 



(72) 



As explained in Sect. 3.1.2 we can also calculate the Compton amplitude as an overlap 



of LCWFs in the symmetric frame of Sect. |2.2j . Using the same approximations as in the 
handbag calculation, Sect. P~T| , and comparing with fl3"0|), (|3~T|) we obtain the analogues 
of the Drell-Yan formula (|66|), (|37l) for our form factors Ry and Ra- With our Gaussian 
ansatz fl4*3|) for the k^-dependence of the LCWFs and the integral ( |5lD the form factors 
Ry and Ra can then be written as 
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+ e 2 d [Ad { v N) (x) + 2 Ad (iV) (x)} + e 2 s 2 As {N) (x)} 



(73) 



Our numerical 



in close analogy to the expressions (^) for the Dirac form factor F\. 
predictions for Ry and Ra are shown in Fig. || 

If one makes the assumption that Ra = Ry as was done in |17| then one obtains the 
suggestive result that the cross section for Compton scattering on the proton is just given 
by the familiar Klein-Nishina expression for Compton scattering on a free quark times the 
square of the form factor Ry(t), which describes the target structure. In our model the 
ratio of Ra and Ry is however rather far from 1 for the values of —t we consider, cf. Fig. ^|. 
From (|73|) one sees that Ra ~ Ry would require all quarks and antiquarks to be completely 

AqW(x) for all N and a, in the range of x 



(AO, 



x 



polarised along the proton spin, i.e. q^ J 
dominating the integrals. For u-quarks this holds indeed if x is close to 1, but not for the 
intermediate x that are important at our values of —t, while for d-quarks the unpolarised 
and polarised quark distributions even have opposite sign. We therefore keep both terms 
R v (t) + R 2 A (t) and R v {t) — R\{t) in (^|); they reflect the fact that the proton target has 
a nontrivial quark spin structure. Using measurements at different values of s and t and a 
Rosenbluth-type separation it will in principle be possible to isolate the new form factors 
|.Ry(t)| and 1^(^)1 from sufficiently accurate experimental data, and to compare them 
with our predictions. 

In Fig. [7| we show our results for the Compton cross section. Given the quality of 
the data, and the small energies and low values of —t and — u at which they are avail- 
able, our predictions compare fairly well with experiment. As a minimum condition for 
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Figure 6: The form factors Ryif) and -R^(t), evaluated with our model for LCWFs. 



our approximations discussed in Sect. fO] to be applicable we only take into account data 
points satisfying —t, —u > 2.5 GeV 2 . Better data and data at larger energies are definitely 
required for a severe check of the new approach and its confrontation with the hard scat- 
tering mechanism. For comparison we also show predictions for the Compton cross section 
at a photon energy of 12 GeV that may be reached at an upgraded JLab facility [53]. At 
such an energy and at cm. scattering angles around 90° the kinematical conditions for 
the approach presented here would be satisfied. Still higher energies, perhaps accessible at 
ELFE |44| , would be even better. 

Dimensional counting ||6| predicts that the Compton cross section scaled by s 6 only 
depends on the ratio t/s or, equivalently, on the scattering angle 9 in the photon-proton 
cm. From Fig. [7| one observes that the soft contributions do not exhibit this counting 
rule behaviour, although they are close to it. s 6 -scaling holds in our approach as long as 
Rv and Ra behave as 1/t 2 . As one can observe from Fig. |] this is approximately true 
for —t in the range from about 5 to 15 GeV 2 . For energies between, say, 3 and 6 GeV in 
the laboratory frame such i-values are only reached in the backward hemisphere. In this 
case the energy dependence of the scaled Compton cross section is hence much milder in 
the backward than in the forward hemisphere (see Fig. [?p. For energies as large as for 
instance 12 to 15 GeV the situation is reversed. The t-values are so large in the backward 
hemisphere that Ry and Ra do not behave as 1/t 2 any more but gradually turn into the 
soft physics asymptotics 1/t 4 . Consequently the scaled Compton cross section exhibits a 
stronger energy dependence in the backward hemisphere than in the forward one. For very 
high energies the soft physics contribution to the large angle Compton cross section scales 
as s~ 10 . We note that Radyushkin's result [17 that all the curves for different energies 
intersect each other at cos 9 = —0.6 does not hold in general. This result may depend on 



specific assumptions made in [17| and holds at best in a rather limited region of energy. 
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Figure 7: The Compton cross section scaled by s 6 
angle in the cm. Data, for —t, —u > 2.5 GeV 2 



versus cos 9, where 6 is the scattering 

151 . Left: 



only, are taken from Ref 
Model predictions obtained from the GRV parametrisation |28j for various photon energies 
in the laboratory frame. Right: Model predictions decomposed into separate Fock state 
contributions at a photon energy of 5 GeV. 



It is an goal of utmost importance to test the energy dependence of the Compton cross 
section experimentally in the relevant kinematical region s, —t, —u ^> m 2 . The present 
data are neither accurate enough nor really satisfy the kinematical requirements. 

The Compton cross section has also been calculated within perturbative QCD P7 | 



and within a diquark model [2y| that combines perturbative elements with additional soft 
physics (correlations in the proton wave function modelled as diquarks). Both models can 
also account for the data although, as we said before, the quality of the data is insufficient 
for a severe test of the models. The diquark model does not lead to the dimensional 
counting behaviour either; it turns out that the energy dependence of the scaled cross 
section in the forward and backward hemisphere predicted by that model is opposite to 
the one of the approach proposed here and shown in Fig. [!]. 



In the leading twist hard scattering calculation of [0 proton distribution amplitudes 
are employed which are strongly concentrated in the end point regions, and thus differ 
drastically from the one determined in Ref. |16| and used here (cf. Eq. (f42|)). For such 
distribution amplitudes the perturbative analysis of Compton scattering, quite like that of 
the nucleon form factor || , may be afflicted by large contributions from the soft end point 
regions, where perturbative QCD is not applicable as we mentioned in the introduction.^ 

19 Rv(t) and RA{t) evaluated from a wave function composed of the distribution amplitude proposed 
in Ref. Q and the Gaussian (|4^) exhibit approximate 1/t 2 scaling behaviour in a much larger t-region 
than found from the distribution amplitude (^)- Also the maximum values of Rv{t) and RA(t) are larger 
by a factor 5 to 8. This parallels the behaviour of the electromagnetic form factors, see Ref. [fill. As 
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We emphasise that in the perturbative approach the dimensional counting rule behaviour 
of the Compton cross section is modified by powers of log s arising from running of the 
the strong coupling constant a s (da/dt oc cef) and from the evolution of the proton wave 
function. These effects have not been taken into account as yet. It remains to be seen how 
much these logs will change the results quoted in Ref. ||47|| . One may also expect that the 
inclusion of transverse momentum effects and Sudakov suppressions in the perturbative 
analysis leads to a similarly strong reduction of the Compton amplitude as was found for 
the proton form factor |12[. In view of this it seems premature to us to claim a success of 
the purely perturbative approach in Compton scattering. 



9 Skewed parton distributions 

In this section we are going to investigate skewed parton distributions (SPDs) || pEB| . 
These distributions are the non-perturbative input for Compton scattering in the deep 
virtual region of small —t but large Q 2 and s. Factorisation of the process into hard and 
soft physics || assures that, like the usual parton distributions, the SPDs are universal 
in the sense that they occur in different hard processes, e.g. in hard meson production. 



As explained in Sect. |3.1.2| we will restrict our investigation of SPDs to the kinematical 
region where they describe how a parton with momentum k is taken out from the proton 
with momentum p and, having undergone a hard scattering, inserted back with momentum 
k + A as a parton inside the scattered proton with momentum p + A (see Fig. 0(b)). 
Due to this restriction we are unable to calculate the full amplitude of the deep virtual 
Compton scattering process, which includes the region < x < ( where we do not have 
a simple representation of SPDs as an overlap of LCWFs. In a restricted kinematical 
region, however, we are able to calculate the process independent SPDs, which are of 
interest in their own. We will also be able to check whether they behave correctly in 
the formal forward limit A = 0, and whether they satisfy bounds coming from positivity 
requirements [|9|, ^(J, [51]]. 

To date essentially nothing is known experimentally about skewed distributions. How- 
ever, various model estimates of the SPDs have been made recently: for instance a bag 



model calculation |52| , a chiral quark-soliton model , and a scalar toy model |3] . A num- 



ber of simple ansatze has also been proposed [54]. In particular the question whether 



there is a strong dependence on the skewedness parameter ( is being debated. 
The spin independent skewed distributions are defined by[^] 

p+ f^Ljxp+z- ( p '|^ a (0) 7 +^( 2 -)|p)=^(x;t)^') 7 + «(p) + 'X-term", (74) 

J Z7T s 

a consequence the Compton cross section does not show approximate s~ 6 scaling behaviour for photon 
energies between, say, 3 and 15 GeV. 

20 For convenience we do not display the link-operator needed to render the definition gauge invariant, 
assuming the use of a light-cone gauge combined with an appropriate choice for the integration path which 
reduces the link-operator to unity. 
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where here and in the following a denotes a quark flavour, antiquarks being explicitly 
labelled by a. The /C-term in ( j74[) goes with the tensor current of the proton and is related 
to proton helicity flip. Like the Pauli form factor F2 and our form factor Ft? we cannot 
evaluate it in our model as explained in Sect. |4.2| . In the definition (|74|) we follow the 
conventions of Radyushkin for nonforward distributions, cf. Eqs. (9.1) and (9.2) of |J. The 
kinematical variables x and ( turn out to be most convenient for calculating the overlap 
of LCWFs. The relation to Ji's original definition of off-forward distributions, where a 
different choice of variables is made, can be found in Ref. ||, Eqs. (9.6) and (9.7), and in 
Ref. [|U Eqs. (24) and (25). 



The matrix element in ([74]) is nonzero in the range — 1 + ( < x < 1 , cf . 0, |25], [18] . 
Re-interpreting a quark with negative momentum fraction as an antiquark with positive 
fraction one finds that — 1 + ( < x < describes the emission and absorption of an 
antiquark, just as ( < x < 1 does for a quark, while in the region < x < ( the proton p 
emits a quark- antiquark pair and is left as a proton with momentum p + A. 

The definition ( ff4| ) reveals the close relationship of SPDs with the usual quark distri- 
butions and with the Dirac form factor. Indeed one finds the reduction formulas 

F ( a =0 (x;t = 0) = q a (x) (75) 

and 

E^T F ( a (x;t)dx = F 1 (t). (76) 

Eq. (^) can be explicitly checked in our results, while we cannot evaluate the moments in 
(|76|), which contain the region < x < ( we do not model here, except in the case ( = 0. 

We now turn to the derivation of an overlap formula for the SPDs. In close analogy to 
the steps that lead to Eq. (|29"D the amplitude for DVCS can be written in terms of proton 
matrix elements as 

a = EKJ 2 f -^== E / T 1 elxp+z ~ M &(°) ^ +1J ^ ^ z ~) \p) 

7 -!+C 2,J\xx'\ A 2 
x [6(C <x<\) u(k', X)H(k', k)u(k, A) 
- 9(0 < x < C) v(-k', -X)H(k', k)u(k, A) 

+ 0(-l + C < x < 0) v(-k', -X)H(k', k)v(-k, -A) 1 , (77) 



with the conventions for spinors given before Eq. (p8|). The different kinematical regions 
mentioned above can easily be recognised. The hard scattering is now approximated as 
collinear, neglecting —t and m 2 compared with Q 2 and setting k = [xp + ,0,0j_], k! = 
[x'p + , 0, OjJ. On the other hand, direct calculation of the overlap diagrams starting from 
the Fock state decomposition of the proton (cf. Sect. |5|) gives the contribution of the region 
C < x < 1 to the amplitude as 

A' = EK,) 2 /'dxEEE [[dxUd^usix-x^^L^ (1-c) 1 ^ 



C N j ft \JXjXj 

x Vlrpi&i, k' ±i ) y N/ 3(xi, ku) u(k', X)H(k', k)u(k, X) , (71 
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where j runs over quarks of flavour a. Note that the label (3 includes a dependence on the 
parton spin A. The arguments of the outgoing wave function x[ and k' ±i are related to Xi 
and k±i by fllCf). From the comparison of ([78D with (ff7[) and the definition (|T4|) we obtain 
the overlap formula for spin-independent SPDs in the region ( < x < 1: 



a (AT), 



C 



x; t) 



CY 



[cLc]jv[d 2 k_L]jv 8(x 



x j, 



iV/9 



EE 

i P 

(79) 

with j again running over all quarks of flavour a. Comparing with fl56|) we see that the 
boundary condition ( [75] ) is correctly implemented in our approach. As for the sum rule 
(ff6|) we find with (|I5D and (|16|) that in the case ( = the overlap expression flT9|) and the 
corresponding contribution from antiquarks reproduce the Drell-Yan formula (|66|) . 

We notice that for ( < x < 1 the r.h.s of ( [79|) has the structure of a scalar product 
in the Hilbert space of wave functions ^Np(xi, k_u). Writing down the Cauchy-Schwarz 
inequality for ([79]) and using the reduction formula ( |75|) we find 



a (AH 



X] t) 
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qi JV) (x) 



Q (iV) 



id) 



for the contribution of each Fock state. Notice that at the points x = ( and x = 1 both sides 
of ([80]) are zero because in the corresponding overlap integrals there are wave functions 
taken at their end points. Summing (|79|) over all Fock states one obtains the analogue of 
(p0|) for the complete distributions; it is precisely the positivity constraint on SPDs derived 
by Pire, Soffer and Teryaev |50|] , which is thus satisfied by the overlap formula ([79]) .^ 

To discuss the emission and reabsorption of an antiquark it is useful to define ^(x; t) 
by the r.h.s. of (|74|) with the field operators replaced with the charge conjugated ones. 
One easily finds the relation J r l ? , (x;t) = —J 7 ®^ — x;t). In the region £ < x < 1 the 
distribution T£{x\ t) describes the emission of an antiquark with momentum fraction x 
and its reabsorption with fraction x' = x — (; along the same lines as above one obtains 
its overlap representation as the r.h.s. of ( |7"9"P with j running over antiquarks instead of 
quarks. One then has of course the analogues of the reduction formula (|75|) and the bound 
(p0|) for jF^ a (x;t) and the usual antiquark distributions. 

Inserting our iV = 3, 4, 5 Fock state wave functions of Sect. |5] in (|79|) we obtain for the 
skewed u and d valence distributions 
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21 Notice that it is satisfied for all t in the physical region t < —( 2 m 2 /(l — £), cf. (fj), with the upper 
bound being i-independent. 
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given at the end of Sect. [| and make use of the relation ( |T0D between Aj_, t and (. The 
exponents n a , m a (N) and the coefficients c^' and are the same as for the valence 
quark distributions discussed in Sect. [| cf. Eq. (|61|) and Tab. |2|, and l g is the number of 



gluons in the corresponding Fock state. For £ = our result (pl|) simplifies to 

1 — x 



X 



(82) 



which is the origin of our simple representations (|68|) and (|73|) of form factors. Finally 
we find that with our wave functions the skewed antiquark and s-quark distributions are 
related with the d valence distribution by the analogues of Eqs. (|59| ) and (|6Q]) . 

In Fig. H| we display our results ( pi]) summed over the N = 3,4,5 Fock states for 
fixed t- values of —0.5 GeV 2 and —1.5 GeV 2 . We remember from the end of Sect. that 
—t > ( 2 m 2 /(l — (). At fixed t this imposes ( < ( max with 



Jt (t - 4m 2 ) + t 

Cma, = 1 ^ • (83) 

We remark at this point that the t-dependence of the SPDs residing in the factor T # does 
not factorise in our approach but mixes with the dependence on x and ( in the exponent of 
([53]); note that the transverse momenta k.' ±i in the overlap formula (|7^) implicitly depend 
on x, ( and Ai through (|16|). A significant dependence on the skewedness parameter ( 
shows up in our results; a fact which is not surprising since ( determines the momentum 
fraction of the active parton in the light-cone wave function of the outgoing nucleon. 

In Fig. |9| we plot the skewed u valence distributions at fixed t and ( as a function of 
x, comparing the contribution from the N = 3 Fock state with the result summed over 
N = 3,4,5. As for the usual parton distributions we see how higher Fock states become 
more and more important as x decreases. We notice that the values of x where this 
happens increase somewhat with (; this can be understood from the fact that at a given 
x the momentum fraction of the parton going back into the proton decreases with (. The 
area under a curve in Fig. ||] gives the u-quark contribution of the regions — 1 + C<^<0 
and ( < x < 1 to the form factor sum rule ([761). We can see that higher Fock states become 
less important as —t increases, in agreement with what we have found for Fi(t) in Sect. [7|. 

For the usual parton distributions we know that both q v (x) and q(x) become singular 
for x — > 0, which cannot be obtained from any finite number of Fock state contributions, 
all of which vanish at x — 0. The question what the situation is for x —>■ ( in skewed 
distributions, when the momentum fraction x' becomes zero while x remains finite, cannot 
be answered in the framework of this paper. We therefore do not claim that our results 
for the contribution of the first tree Fock states describe the full distribution as x comes 
close to (. 

The definition of spin dependent SPDs is obtained from ( |74]) by the replacements 7 + — > 
7 + 7 5 and — > </ °; for antiquarks one has Q^[x\ t) = — x;t). The analogue of the 
/C-term goes now with the pseudoscalar current of the proton and is again related to proton 
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Figure 8: Skewed parton distributions for u and d valence quarks in the region £ < x < 1, 
obtained from the N = 3,4, 5 Fock states (P3 = 0.17, P4 = P5 = 0.1). The values of ( max 
are 0.52 at t = -0.5 GeV 2 and 0.71 at t = -1.5 GeV 2 . 
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Figure 9: Skewed valence u-quark parton distributions in the region ( < x < 1 as a function 
of x at fixed ( and t. Full lines show the sum of contributions from the N = 3,4, 5 Fock 
states and dashed ones the contribution from the N = 3 Fock state alone. The vertical 
lines in the plots indicate the value of x where the N = 3 contribution is 80% of the one 
summed over N = 3,4, 5. 
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spin flip. From the appropriate overlap formulae we find the spin dependent skewed valence 
distributions 



g ( a{N \x-,t) - g- a{N \x;t) = A6f ) p N (i - cr^-'' r N (x,c, -t(i-0-C 



2 2 y 
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34) 



where the coefficients Afe^, Ac^, and Ac?W are the same as the ones for the spin 
dependent valence distributions listed in Tab. |3|. Evidently the spin dependent skewed 
distributions reduce correctly to the usual ones in the limit £ — > and £ — ► 0. 



10 Summary 

In the present paper we have linked ordinary and skewed parton distributions to soft 
overlap contributions to elastic form factors and to large angle Compton scattering using 
nucleon light-cone wave functions. 

We have investigated how and under which conditions overlap contributions to exclusive 
processes can be expressed in terms of LCWFs. For large angle Compton scattering, at 
large values of the Mandelstam invariants s, —t and —u, we can calculate the soft overlap 
contribution using its factorisation into handbag diagrams, i.e. into soft parton emission 
and reabsorption by the nucleon and a hard parton-photon scattering. In the case of 
deeply virtual Compton scattering, with large Q 2 and s but small —t, where we cannot 
express the amplitude as an overlap of soft LCWFs, we have calculated the skewed parton 
distributions in the large-x region. 

For the LCWF of the three-quark nucleon Fock state we have taken over the parametri- 



sation of [16|, which involves only two parameters adjusted to data. For the Fock states 
with an additional gluon or quark- ant iquark pair we have taken a very simple ansatz, in- 
troducing only two more parameters, which are fitted to the gluon and sea quark parton 



distributions from the GRV analysis | [2B|| . The values of all four parameters come out in 
a range compatible with their physical meaning of Fock state probabilities or a transverse 
size parameter. In the overlap contributions to Compton scattering and the form factor we 
also estimate the net effect of all higher Fock states, using as input the difference between 
the GRV parton distributions and those calculated from the three lowest Fock states only. 

The phenomenology we can do with our ansatz is very rich: we reproduce well the 
unpolarised and polarised parton distributions down to x around 0.5, as well as the data 
for the nucleon Dirac form factors and for real Compton scattering at large cm. angles. 
The inclusion of higher Fock states in the soft overlap contributions confirms that as —t 
increases the lowest Fock states become increasingly dominant and gives an impression of 
the accuracy one can hope for by only taking into account the three quark state. The 



LCWF of UToTI was constructed so as to saturate the elastic form factor data. The fact that 



with the same wave function one obtains a reasonable description of Compton scattering 
supports the hypothesis that there is no sizeable perturbative contribution to either process 
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in the range of momentum transfers where data exist; soft physics seems to dominate as 
was occasionally suggested in the literature [|7|, §, |i~0|. 



We stress that from the apparent agreement of exclusive data with dimensional counting 
rules the dominance of perturbative QCD contributions cannot be deduced. Soft physics, 
as for instance the overlap-type contributions which we propose, provides broad maxima 
in scaled observables such as t 2 Fi{t) and the scaled Compton cross section s 6 da/dt, and 
thus mimics dimensional counting rule behaviour in a certain range of t. 

Compared with the elastic form factors large angle Compton scattering has a second 
independent kinematical variable and thus provides an additional handle to experimen- 
tally test how well dimensional counting rules are satisfied. We further suggest that the 
imaginary part of the scattering amplitude, which is accessible in virtual Compton scat- 
tering with a polarised lepton beam, offers a sensible tool to investigate which dynamical 
mechanism is at work: in the handbag mechanism imaginary parts are only generated 
through loop corrections to the photon-parton subprocess, whereas in the hard scattering 
mechanism real and imaginary parts generically are of the same order of magnitude. Spin 
observables may also be sensitive probes of the underlying physics, given the particular 
helicity structure of the photon-parton scattering in the handbag diagrams. In any case 
we see a strong motivation to have further and more accurate Compton data at sufficiently 
high values of energy and momentum transfer. 
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